On the generation of density perturbations at the end of inflation 
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Recently a mechanism was proposed whereby the primordial density perturbations are generated 
at the end of inflation. We continue the analysis of the proposed model of this mechanism and 
calculate the maximum extent to which the density perturbations produced via this model can 
dominate over those of the standard inflationary paradigm. In addition, we provide a straightforward 
variation of this model which allows for greater amplification of the density perturbations. Finally, 
we show that a variation in the implementation of the original model results in significant non- 
gaussianities in the resulting spectrum of density perturbations. The level of non-gaussianities can 
be made to saturate the current observational bound. 
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I. INTRODUCTION 

Measurements of the cosmic microwave background ra- 
diation [l| have revealed a highly uniform background 
with relatively scale-free super-horizon perturbations on 
the order of a few parts in 10^. A possible source of 
these perturbations is found in the quantum fluctuations 
of one or more light scalar fields during an early epoch of 
infiation [1,3. This is because in the (quasi-) de Sitter 
space of infiationary expansion, a quantum fiuctuation in 
a scalar field evolves according to classical equations of 
motion after its wavelength exceeds the rapidly decreas- 
ing Hubble length 0] . A fiuctuating scalar field may then 
be converted into an energy density perturbation via a 
variety of proposed mechanisms. 

For example in the standard inflationary paradigm fsj , 
inflation is driven by the potential energy of a single 
slowly-rolling scalar field. In this case the infiaton can 
be viewed as a unique clock parameterizing the evolution 
of the early universe. Therefore fiuctuations in the infia- 
ton translate into fiuctuations in the duration of inflation. 
Since the energy density of the universe redshifts more 
rapidly after inflation than during inflation, this results in 
energy density fluctuations on surfaces of constant scale 
factor after inflation. 

More recently, it was proposed that energy density per- 
turbations could result from fluctuations in light scalar 
flelds that do not contribute significantly toward the in- 
flationary dynamics. For example, in the "curvaton" sce- 
nario (see also [^) a light scalar field dubbed the cur- 
vaton receives fiuctuations. After infiation, the curvaton 
evolves as a massive fluid and therefore redshifts more 
slowly than the radiative products of reheating. If the 
curvaton eventually dominates the energy density of the 
universe and then decays, density perturbations result 
because the duration of curvaton domination depends on 
the fluctuating curvaton. The inhomogeneous reheating 
scenario [7] also achieves density perturbations by vary- 
ing the duration that a massive fluid dominates the en- 
ergy density of the universe. In this case the duration 



of domination is modulated via a decay width that de- 
pends on some fiuctuating light scalar field. Other ways 
to modulate the duration of of massive fluid's domination 
have been proposed in Q . These ideas have recently been 
reviewed in [9|. 

Another major variant on these proposals was recently 
offered in ^lOj (note however that a similar idea was pre- 
viously used in [01). Whereas in the standard picture 
the duration of inflation is influenced by the fiuctuations 
in the infiaton as relevant modes leave the horizon, in this 
picture the duration of inflation is influenced by fiuctua- 
tions in the infiaton when inflation ends. This can occur, 
for example, if inflation ends at an inflaton value that 
depends upon some other field that receives fiuctuations 
during inflation. In a specific model was used to 
demonstrate that this scenario can lead to significant am- 
plification of the density perturbations that result from 
slow-roll infiation. 

Here we continue the analysis of that model and calcu- 
late the maximum extent to which the resulting density 
perturbations can dominate over those of the standard 
infiationary paradigm. We also explore the sensitivity 
of this result to the tuning of model parameters. In ad- 
dition, we provide a straightforward variation upon this 
model which allows for greater amplification of the den- 
sity perturbations. Finally, we show that a slight varia- 
tion in the implementation of the original model allows 
for significant non-gaussianities in the spectrum of den- 
sity perturbations. The level of non-gaussianities can be 
made to saturate the current observational bound. 

This paper is organized as follows. In Section [H] we 
summarize the mechanism of [lo| and establish a conve- 
nient formalism and notation. Section IIIII describes the 
specific model introduced in [l^ , and in Section IIVI we 
analyze this model in greater detail. ScctionlVlintroduces 
two variations upon this model and studies their conse- 
quences. Discussion and concluding remarks are provided 
in Section IVll 
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We parameterize metric perturbations using the cur- 
vature perturbation on surfaces of constant scale factor, 



otherwise known as the Bardeen variable C 
corresponds to writing the interval 



This 



-dt^ + a'^{t)e^'^^''''^-f,j{t,x)dx'dx^ , (1) 



where a and C are chosen to give jij unit determinant. 
In this work we ignore tensor perturbations, which allows 
us to write 7^ (t, x) = dij since scalar field fluctuations 
carry no anisotropic stress (for a review of cosmological 
perturbation theory see p^). 

The proposal of is most easily studied within the 
so-called SN formalism of (see [l5| for recent ex- 
tensions). This formalism reveals that on super- horizon 
scales the Bardeen variable C can be interpreted as the 
perturbation in the number of e-folds of local expansion 
between an "initial" spatially-flat hypersurface (at, say, 
time t = to) and the comoving hypersurface on which C 
is evaluated M. 



C(i,x) = iV(t,x) - In 
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On super- horizon scales, spatial gradients can be ignored 
and N{t, x) can be computed by treating large, widely 
separated regions as separate FRW universes, but with 
different initial field values [l4j • 

We limit our attention to the case where inflation is 
driven by a single inflaton. Then fluctuations in N come 
only from the value of the inflaton when its wavelength 
becomes larger than the Hubble length, denoted (j)k, and 
the value of the inflaton at the end of inflation, denoted 
(pc- Thus the number of e-folds of universal expansion is 
iV((/)fc,(/)o), and 
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That the power spectrum for C is observed to be predom- 
inantly Gaussian allows us to neglect higher order terms 
that might appear in the expansion of Eq. ([3]) . 

In the standard picture, inflation ends when the slow- 
roll conditions are violated. For single field inflation this 
happens at a unique value 4>e and therefore the second 
term in Eq. ([3]) is zero. For inflation to end at varying 
values of (jje requires the addition of at least one other 
field. For example, with an additional field a it is possible 
that (/'c(o') depends upon a and inflation ends within a 
range of 0o given by 
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(4) 



If the masses of the scalar fields are much less than the 
Hubble rate the fields acquire a constant power spectrum 
of fluctuations Vstjik = 'Psak = (-fffe/27r)^, where Hk is the 
Hubble rate at the time the mode k exits the horizon [l^ . 
Thus the power spectrum for C. is 



The density perturbations generated at the end of infla- 
tion will dominate over those produced from the standard 
picture when 



dN . 



da 



> 



ON 



(6) 



Note that dN/d(j)k 7^ dN/d(j)e since the former probes 
the dependence of TV on </) deep in the inflationary epoch 
while the latter probes the dependence of on ^ near 
the end of inflation. 



III. THE SPECIFIC MODEL 

The above ideas were demonstrated in [lo| using a spe- 
cific model described by the potential 
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+ \\a(T^X^ + V,{a). 



(7) 



As was noted in ^lO], this is the original hybrid inflation 
model 17] but with interactions involving an additional 
light scalar field a added with the last two terms. 

Hybrid inflation assumes the initial conditions of 
chaotic inflation [13, [3- However, during the early 
stages of hybrid inflation the field x is pushed to zero 
much faster than (and likewise faster than a in the 
above model). Then (p becomes a slowly- rolling inflaton 
withjnflation assisted by the vacuum energy m^/Ag. As 
in [101] we assume that the self interaction of a repre- 
sented by the term Va- does not contribute significantly 
toward the infiationary dynamics. 

The field x is pinned to the origin until the curvature 
in its potential, d'^V/dx^i becomes negative. In |lo| a 
scenario is described in which inflation ends abruptly as 
X rolls away from the origin. This happens when 
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(8) 



Therefore if a receives fluctuations 5a^ infiation ends at 
field values (p^, that vary according to 



'-5a. 



(9) 



Referring back to Eq. ([5]), we see the density perturba- 
tions resulting from fiuctuations 5(j)e dominate over those 
resulting from 5(j)k when 
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In the second equation we have used that the first slow- 
roll parameter can be written 
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An important consequence of this mechanism is that 
the spectral tilt is given by the tilt in the spectrum of fluc- 
tuations Sa, as opposed to ^0. This gives a tilt which is 
independent of the second slow-roll parameter 77 [lol . [l^ . 
Therefore rj is not directly constrained by observation. 
In addition, since e ~ 2He (2e — rj) it appears as if e 
may decrease significantly during the course of inflation 
if ?7 ^ e. Thus it was suggested in [10] that the condition 
of Eq. pO)) is easily satisfied. 

However, one might note that R is proportional to 
Efe which according to this mechanism is constrained by 
observational bounds on the spectral tilt. These give 
Efe ^ 0.02 [l| in the absence of cancellations [13] ■ In ad- 
dition, we do not expect ec and 0c to be independent of 
each other. Thus it is worthwhile to investigate Eq. (fTUl) 
in greater detail in order to determine precisely what lim- 
its the extent to which the density perturbations of this 
model can dominate over those produced via the stan- 
dard inflationary paradigm. 
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where mpi denotes the reduced Planck mass. The first 
constraint of Eqs. (fT4|) allows us to write 
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Plugging these values into Eq. pH)) gives 
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IV. A MORE DETAILED ANALYSIS 

We now consider the model of Eq. ([7]) in greater de- 
tail in order to clarify the constraints on R as given by 
Eq. pn]) . Consistent with the analysis of [10] we ignore 
the contribution of Va toward the vacuum energy. As 
described above, the x field rolls away from the origin 
when 6 takes the value 



A. 



(12) 



To explain the second approximation, first note that we 
require Aq-cTc ^ ^'^ order to avoid x remaining trapped 
at the origin after reaches zero. Meanwhile, the precise 
value of CTo is a stochastic variable constrained by con- 
ditions independent of \a and m-^. Thus to ensure that 
the desired dynamics are typical of this model requires 
we set m^/Acr ^ cr^ for typical values of CTo. 

For inflation to last until (j) — (jy^ but end abruptly when 
X rolls away from the origin requires that the term ra^ j 4g 
dominate the energy density at the end of inflation and 
that be much greater than the Hubble rate at the 
end of inflation.^ These conditions respectively give the 



The quantity (j)k is the value of the inflaton when the 
mode k exits the particle horizon. For cosmological scales 
of current interest this happens about 60 e-folds of infla- 
tion prior to the end of inflation. 

In order to find the maximum value for R, we find the 
value of (j)k that maximizes R and set the parameters g, 
m^, and such that the mode k exits the horizon 60 
c-folds prior to inflation. This is equivalent to finding 
the balance between the vacuum energy ni'^/Ag and the 
potential energy m?0^ that maximizes R. The result is 
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with the constraint that the mode k leaves the horizon 
at Nk~ 60, with 
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In the last expression we have used the constraints of 
Eqs. (fT4)) to identify the most significant term. Putting 
this all together we find the maximum value of R to be 



The other possibility is to set (pa = n^x/\/ ^<t> *° '-'^ greater than 
the value of </> for which inflation ends when the potential of (f> 
dominates. Since in this case inflation ends when </> is of the 
order of the Planck mass Il8l . this requires a Planck scale 
unless S> 1. Nevertheless, an analysis similar to that which 
follows gives 
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where 0j, is 4> evaluated 60 e-folds before the end of inflation. It 
will be seen that this is no better than the result obtained in the 
scenario analyzed in greater detail in what follows. 
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Before proceeding to study Eq. ([^0]) we should check 
that producing the correct power spectrum normalization 
does not introduce any constraints that conflict with our 
present assumptions. Applying Eq. ([5]) to the scenario 
considered above we find that when density perturbations 
generated at the end of inflation dominate we have 
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To match observation we require Vq to be very small, 
Vc^ « (5x10-^)2 According to the second of Eqs. ^ 
the first term in parentheses is already constrained to be 
much less than order unity and in fact can be set as 
small as necessary to match observation. In addition we 
expect / rripj to be very small. Thus to set Vc, to match 
observation does not introduce any constraints in conflict 
with those of the above analysis. 

We have written R in the form of Eq. (pO| in order to 
emphasize the maximum extent to which the perturba- 
tions produced at the end of inflation may dominate over 
those produced as cosmological scales exit the horizon. 
The first term is weakly constrained by the energy scale 
of inflation and is here taken to be Nu « 60. The first 
term in parentheses is constrained to be significantly less 
than unity as explained in the discussion below Eq. (|12p . 
Finally, the argument of the logarithm must be much 
greater than unity in order to satisfy the first constraint 
of Eqs. (HH). Therefore the last term is less than unity. 
Thus for appropriately tuned parameters we might ex- 
pect the product of these three factors to be a couple 
orders of magnitude below unity. 

However, Eq. (|^ also contains a factor rn^Jal. In 
(quasi-) de Sitter space a scalar field such as a evolves 
both according to its classical equation of motion and 
due to quantum fluctuations as modes leave the parti- 
cle horizon. The net effect of this evolution is that the 
correlation function (cr^) migrates toward a fixed value 
depending upon Va and the Hubble rate H [l^. For ex- 
ample, if inflation lasts long enough and if Va = ^m^a^ 
then (ct^) ^ H'^ /nri^ ^M- Taking a typical value of CTo to 
be tTo ~ gives 

The dynamics described above require <^ <C H^, 
so the first term in Eq. must be at least a few orders 
of magnitude below unity. However, the second term can 
be very large. Current observation gives rri^^/H'^ > 10^ 
[20| . which is more than sufficient to compensate for all 
the small factors in R if parameters are tuned appropri- 
ately. Reducing the scale of inflation allows for greater 
values of R. Of course, our Hubble volume could also 
be a region of atypically small CTc- Finally, we may re- 
duce (Tc to an arbitrarily small scale by considering a to 
be a p seudo-Nambu-Goldstone boson (see for example 
BUI). 2 Thus in a number of circumstances we expect 



^ An interesting scenario involves a pseudo-Nambu-Goldstone bo- 
son that ranges over a scale ~ ^x- Then Ao- may be a coupling 
of order unity and we obtain 

R~Nk^, (23) 

where we have dropped the logarithm and other factors of order 
unity. Clearly R is much greater than unity in this case. 



the level of density perturbations generated at the end of 
inflation to be significantly larger than those produced 
when cosmological scales exit the horizon. 

In Eq. (|20p the parameters g, m^, and are tuned 
such that scales of cosmological interest leave the hori- 
zon when (j)k is given by Eq. (|18[) . To study the result 
of relaxing this tuning requires to invert Eq. ()19|) to ob- 
tain 4>k{Nk) and insert the result into Eq. (I17|) to obtain 
R{Nk)- This allows for the remaining parameters in R 
to be freely varied while R retains its original meaning; 
that is, that R compares density perturbations produced 
at the end of inflation to those produced when cosmolog- 
ical scales exit the particle horizon. 

Note that the first two terms in the brackets of Eq. (fTg]) 
always dominate over the third term and that they are 
comparable to each other when (j)k is given by Eq. (|18|) . 
Remember that (j>k is defined as the value of cj) ^fe ~ 60 
e-folds before the end of inflation. Decreasing slows 
the evolution of which therefore decreases (j)k- In this 
case the first term in brackets becomes more important 
and inverting Eq. gives 

^l^c^le.p(8N,'-^y (24) 

For less than in Eq. (fT8|) the important functional 
dependence of R is R{4>k) oc 0^. Thus we see that R 
decreases exponentially when the ratio gm^^m^^/m^ is 
decreased from its optimal value. 

On the other hand, increasing quickens the evo- 
lution of (f) and therefore increases 4>k- In this case the 
second term in the brackets of Eq. becomes more 
important. Inverting in this case gives (/)| « Nkiri^^ 
which is relatively independent of the model parameters. 
Therefore the magnitude of R changes predominantly 
through its dependence upon m'^(j>f. in the denominator 
of Eq. (I17|) . Thus we see that significantly increasing 
the ratio gm'^mp^/m'^ from its optimal value results in a 
roughly proportional decrease in the in size of R. 

V. GENERALIZING THE MODEL 
A. Varying the Potential for (j> 

According to Eq. ([6]) , the generation of density pertur- 
bations at the end of inflation is most effective when the 
slow-roll parameter near the end of inflation is much less 
than when cosmological scales of interest exit the particle 
horizon. Thus to enhance the resulting perturbations we 
desire a potential for (p that decreases more steeply for 
deep in the inflationary epoch and decreases more gently 
for (f) near the end of inflation. This can be accomplished 
by replacing the ^m'^cj)'^ term in Eq. ([7]) with jXcj)'^. 

We analyze this scenario in exact analogy to the anal- 
ysis in Section IIVI Again we assume that the vacuum 
energy at the end of inflation is dominated by th e /Ag 
term. The calculations proceed just as in Section HvT and 
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in the end we find 



(25) 



As described in Section IIVI the term in the first set of 
parentheses is expected to be significantly less than unity. 
However the last term is expected to be greater than 
unity in order that x becomes pinned to the origin in the 
early stages of inflation. In addition, the numerical pre- 
factor IQSNl - 10^ for Nu ~ 60. The discussion about 
the factor ra^^/al near the end of Section HVl applies here. 
Thus we see the maximum level of density perturbations 
produced at the end of inflation can be greatly amplified 
by simply introducing a 0^ potential. 



B. Relaxing the Constraint on \a 



where C,g symbolizes a variable with a gaussian spectrum 
[2^ . In this case is the curvature perturbation pro- 
duced via the standard inflationary paradigm, 



(28) 



Here we have assumed for simplicity that scales of cosmo- 
logical interest exit the horizon after has rolled to zero. 
Thus we treat the relevant stages of inflation as driven 
by the a field, assisted by the vacuum energy m^/Ag, 
with (j) being a heavier field fluctuating about the origin. 
We neglect the non-gaussian component of C, produced 
when relevant scales exit the horizon as this is in general 
relatively small [l^. For comparison to observation it 
does not matter whether the non-gaussian fluctuations 
are sourced by the same field as the gaussian fiuctua- 
tions [111. Therefore we can combine Eq. ([3]), Eq. ((271) 
and Eq. ([28]) to obtain 



In Section Hm we followed the assumption presented in 
[l0| that near the end of infiation \c<J^ < rn^ so that 
infiation ends before (f> reaches zero. It is interesting to 
explore the consequences of lifting this assumption. If we 
do not demand that X^a^ < rn^, then for large enough 
A(j or a the x field is still pinned to the origin when (j) 
reaches zero. We then expect infiation to continue with 
the a field rolling down its potential until the condition 
of Eq. ([5]) is met. At this point x I'oHs away from the 
origin and abruptly initiates the end of inflation. 

Nevertheless, if <^ H during inflation, then 
will retain a power spectrum of fluctuations given by 
7^50^ = {Hk/2-KY. Thus the fields a and (j) have essen- 
tially changed places, with a playing the part of the infia- 
ton and (p the fluctuating fleld. However, in this scenario 
(j) has no homogeneous component and 
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The density perturbations produced at the end of infla- 
tion are now entirely non-gaussian! 

It should be noted here that the original model de- 
scribed in Section IIIII also results in non-gaussianities. 
As described in in that model the non-gaussianities 
result from higher order terms in the expansion of SN 
in Eq. ([3]). The situation we consider here is different in 
that the non-gaussianities arise at leading order in the 
expansion of SN. 

The level of non-gaussianity can be parameterized with 
a term /nl defined according to the equation 



C — Cg — 7/NLCe 



(27) 



^ This subsection fails to account for the lack of correlation be- 
tween a and fluctuations when computing the non-gaussianity 
parameter /nl- The issue is discussed in the erratum. 



/NL 



5 (dN 
3 \dae 



3 Act (To da^ V dak 



(29) 



If we take Va — ^m^a'^, it is straightforward to calculate 



the largest possible /nl by translating the arguments of 



Section [IV 
and that 



In particular, we note that tTe 
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Likewise, minimizing the factor [dN/dok) ^ gives Ok — 
rr? I ^flg-nif, and 
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Finally, putting all this together we find 
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It must be emphasized this is an upper limit; smaller 
/nl are easily achieved by choosing parameters that do 
not minimize {dN / dak)~^ ■ Note also that in this imple- 
mentation of the model the only constraint on A^ is that 
Xtj.ScjP' ~ X^{Hk/\/2Y <C TO^. This constraint assures 
that X always rolls away from the origin before a reaches 
zero, and gives 
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Referring to second constraint of Eqs. (|14p . we see the 
term in parenthesis is already constrained to be much 
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less than unity. Therefore this mechanism permits a non- 
gaussian component to the density perturbations all the 
way through the observational limit of /nl ^ 135 [13] ■ 

For the above calculation to be appropriate requires 
that the Cg of Eq. (|28|) actually be the dominant con- 
tribution to the curvature perturbation. For this to 
be the case its power spectrum must have a magnitude 
Pc,s ~ (5 X 10~^)^ in order to match observation l]. 
From Eq. ([5]) and Eq. (|3T|) we have 
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where we have used that in this scenario the first slow-roll 
parameter as the mode k exits the horizon is 
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We have written the power spectrum in this way in order 
to employ the observational constraints on the spectral 
tilt and on the level of gravity waves. 

Since the primary, gaussian density perturbations are 
now sourced during inflation, the spectral tilt is n ~ 1 — 
2r] — 6efe w 2efe ^ 0.04 [1], where we have used that the 
second slow-roll parameter in the scenario we are consid- 
ering is given by 



V do 
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Meanwhile, it is observed that Hl/rn^^ < 10"^ In 
this model we can decrease the spectral tilt independently 
of H^/rrip^ by simply decreasing m^-. Therefore we sim- 
ply require to set H^/iri^^ ^ 10~^ to match observation. 



Here 
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where we have expressed the result in terms of /nl to 
clarify which are the remaining free parameters. Ac- 
cording to the second constraint in Eq. , the term in 
parentheses is already constrained to be much less than 
unity. In fact, it can be set as small as necessary to satisfy 
observation. Moreover, since the three terms in Eq. ((37|) 
depend upon three independent parameters (A,/,, g, and 
TO^), we have considerable freedom in exactly how we 
satisfy the observational bound. Thus we conclude that 
this model allows for significant non-gaussianities for a 
range of model parameters. 



VI. DISCUSSION AND CONCLUSIONS 

In this work we continue the analysis of the model pro- 
posed in [10] to generate density perturbations at the end 
of infiation. We confirm that these density perturbations 



can easily dominate over those produced via the standard 
inflationary paradigm, and explore the sensitivity of this 
result to the tuning of model parameters. In addition, we 
provide a straightforward variation of this model which 
allows for even greater amplification of the density per- 
turbations. 

It is worthwhile to consider how general is this analy- 
sis. According to Eq. ^ the production of density per- 
turbations at the end of inflation is most effective when 
the slow-roll parameter near the end of inflation is much 
less than that when cosmological scales of interest exit 
the particle horizon. Since inflation can only end with 
the slow-roll parameter rising to unity, this suggests the 
mechanism is most effective only when the inflationary 
potential contains large derivatives near the point where 
it dips toward its minimum. These large derivatives are 
most naturally accomplished by inserting a second field 
direction for the vacuum energy to fall to zero and ini- 
tiate a reheating phase. This is precisely the scenario 
implemented in hybrid inflation and generalized in the 
models considered here. 

We also study a variation in the implementation of 
the model proposed in that results in modifled in- 
flationary dynamics. We show that this case results in a 
spectrum of density perturbations with significant non- 
gaussianities for a range of model parameters. In partic- 
ular it is shown that these non-gaussianities are capable 
of saturating the current observational bound. 



VII. ERRATUM 

In Eq. (j27p we define the non-gaussianity parameter 
/nl and in Eq. (P^)) we relate it to model parameters. 
This definition and the subsequent calculation is appro- 
priate when the gaussian and non-gaussian contributions 
to C come from fiuctuations in the same field, or when 
they are correlated as such, but when they are uncor- 
related, a different approach is more appropriate [25j . 
What is defined as /nl above corresponds to a quantity 



/nl = -o'^C^c 



(38) 



where r<^ is the non-gaussian fraction of Ref . [25| . As ex- 
plained in Ref. 25], the non-gaussianity parameter that 
is relevant to observation is given by 



/nl ^ -rlV^ 



-1/2 



(39) 



where some factors of order unity and a logarithmic de- 
pendence on the scale of observation k have been sup- 
pressed. Combining with Eq. ([55)1 gives 



/i 



NL 
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(40) 



This result is suppressed by a factor "P^^ « (5 x 10^^)^, 
but includes the quantity /nl of Eq. ((5^ cubed. There- 
fore, whereas it is argued in Section IV Bl that the right- 
hand side of Eq. ([5^ need only approach 0(100) to reach 
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the observational limit, we here find it must approach 
0(10'*) to reach this limit. Still, this appears to be easy 
to achieve in this model. Note also that these results as- 
sume the vacuum expectation value of the light field cf) 
is precisely zero (no homogeneous component), as is the 
correlation between and cr. Violation of these assump- 
tions (due for instance to averaging over just one Hubble 
volume) may lead to increased levels of non-gaussianity 
for the same parameter values. 

The author is grateful to Christian Byrnes for explain- 
ing these issues to him. 
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